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Abstract. We consider the quantum site percolation model on graphs with 
an amenable group action. It consists of a random family of Hamiltonians. 
Basic spectral properties of these operators are derived: non-randomness of 
the spectrum and its components, existence of an self-averaging integrated 
density of states and an associated trace-formula. 



1. Introduction: The Quantum percolation model 

The quantum percolation model (QPM) consist of two building blocks which are 
both well studied in physics of disordered media. 

Let us first introduce the classical site percolation model. It is used to model 
the flow of liquid through porous media, the spreading of forest fires or of diseases 
etc. Consider the graph Z d where two vertices are connected by an edge if their 
Euclidean distance equals one. Equip each vertex v G Z d with a random variable 
q(v) taking the value with probability p and "oo" with probability 1—p and being 
independent of the random variables at all other vertices. For each configuration of 
the randomness ui G 57 let V(u>) :— {v G Z d | q(v) = 0}. The percolation problem 
consists in studying the properties of connected components — called clusters - 
of V(uj). Typical questions are: With what probability do infinite clusters exist? 
What is the average vertex number or diameter of a cluster? What is the probability 
that 0, v G 1i d are in the same cluster, etc.? One of the central results of percolation 
theory is the existence of a critical probability p Cl such that for p > p c (respectively 
for p < p c ) an infinite cluster exists (respectively does not exist) almost surely. See 
e.g. [11, 10] and the literature cited there. 

Random lattice Hamiltonians are used to describe the motion of waves in disor- 
dered media. Each of them is a family of operators on i 2 (Z d ) indexed by elements 
of a probability space. The family obeys an equivariance relation with respect 
to the action of a group. More precisely, the group acts in a consistent way on 
l 2 (Tj d ) by translations and on the probability space by ergodic, measure preserving 
transformations. 

The spectral features of these random operators allow one to draw conclusions 
about the transport and wave-spreading properties of the modelled medium. Mono- 
graph expositions of this topic can be found in [4, 16, 19]. 
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Let us define the simplest QPM: Let 
(1) (A ul f){v) = fM for aU v, 10 G V(w) and &\1 f E I 2 (V (uj)) 

dist(v,w) — 1 

be the adjacency operator of introduced above. More precisely, is the 

adjacency operator of the induced sub-graph of Z d with vertex set V(u>). Here 
"dist" denotes the distance function on this graph. 

At this point let us explain why we chose oo as one of the values the random 
variable q(v) takes. The adjacency operator on Z d corresponds (up to an energy 
shift) to the kinetic energy part of a quantum Hamiltonian on the lattice. In this 
picture q corresponds to the potential energy. In the quantum percolation model 
it vanishes on some sites, on others it is infinitely high, i.e. forms an impenetrable 
barrier for the quantum wave function. 

The interesting feature of the QPM is that it defines a Laplacian on random 
geometry. More precisely, its domain of definition l 2 (V(u)) varies with lo. This 
is the main difference to the random lattice operators considered in [4, 16]. After 
an extension of the notion of random lattice Hamiltonians the QPM fits in this 
framework. In our approach we rely on methods from [17, 14, 13], developed there 
to study operators on manifolds. 

The QPM was first studied in [8, 7]. There it was considered as a quantum 
mechanical model for electron-propagation in binary alloys where only one of the 
two types of atoms participates in the spreading of the wavepacket. The model 
attracted special attention because of the existence of molecular states, i.e. eigen- 
vectors supported on finite regions of the infinite cluster, sec [12, 5]. The last cited 
reference is the motivation of the present paper and our results can be seen as a 
mathematically rigorous version of some arguments in [5] and their extension to 
more general graphs. 

The integrated density of states (IDS) of a Hamiltonian is the number of eigen- 
values per unit volume below a certain energy value. Thanks to the stationarity 
and ergodicity assumptions it is well defined for random Hamiltonians. The IDS of 
a random Hamiltonian captures its global spectral features and its understanding 
is the prerequisite of the study of finer spectral properties. In the present work we 
analyze this quantity and provide therewith a basis for a further study of the QPM, 
cf. Section 4 and [20]. 

The next section states the results of this note, Section 2 is devoted to their 
proofs and the last section concludes with a discussion of further research topics. 

2. Results: Spectral properties of finite range hopping operators 

To describe the geometric setting we are working in precisely, let us recall basic 
notions from graph theory and fix the notation along the way. A graph G = (V, E) is 
given by a set of vertices V = V(G) and a set of edges E — E{G) C (V x V\{(v, v) \ 
v G V}) / ~. Here ~ denotes the relation (v,w) ~ (w,v). If e = (v,w) G E, we 
call v,w G V nearest neighbours and endpoints of the edge e. By our definition a 
graph is simple: it contains neither multiple edges nor self-loops, i.e. edges joining 
a vertex to itself. 

A path (of length n) in G is an alternating sequence of vertices and edges 
{vq, ei, Vi, . . . e n , v n } such that ej — (vj-i,Vj) for all j = l,...,n. If there is a 
path between two vertices v and w they are called (path) connected. This relation 
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partitions the graph into (path connected) components. If a component contains a 
infinite number of distinct vertices we call it an infinite component. The distance 
between two vertices v, w E V is defined by 

(2) dist G (u, w) := dist(u, w) 

:= minjlength of p \ p is a path connecting v and w} 

Note that the distance between v and w in a sub-graph of G may be larger than 
their distance in the original graph G. The vertex degree deg(w) of a vertex v E V 
equals the number of edges e G E, such that v is an endpoint of e. 

Let G and G' be graphs. A map <p: G — > G' is called a graph-map or graph- 
homomorphism, if </>: V(G) — > V(G'), (f>: E(G) — > E(G') and if for any e = (u, w) G 
E(G), the image (/>(e) equals 4>{w)). A graph-map <f>: G G which has an 

inverse graph-map is called a graph- automorphism or simply automorphism of G. 

Let r be a group of graph- automorphism acting on a graph X. It induces 
a projection map proj: A — > A/T. We assume that the quotient is a finite 
graph. This implies in particular that the degree of the vertices in V is uni- 
formly bounded. We denote the smallest upper bound by deg + . Chose a ver- 
tex [v] E V(X/T) and a representative v E [v] C V(X). Starting form v, lift 
pathwise the vertices and edges of X/T to obtain a connected set of vertices 
and edges T C X, such that proj|^: T — > A/r is a bijective map. The set 
f := {u £ ^(A) | v is an endpoint of an edge in J 7 } is a graph, which we 
call fundamental domain. Note that proj|^: T — > A/r is a graph-map, which is 
bijective on the set of edges, but not on the set of vertices. 

We construct a probability space (Q,A,F) associated to percolation on A. Let 
^ = x^gylO, 00} be equipped with the cr-algebra A generated by finite dimensional 
cylinders sets. Denote by P a probability measure on fl and assume that the 
measurable shift transformations 

t 7 : Q — > fl, (t 7 uj) v = w 7 -i„ 

are measure preserving. Moreover, let the family t 7 ,7 E T act ergodically on fl. 
By the definition of r 7 ,7 E T the stochastic field q: ft x V — > {0, 00} given by 
g(o;, u) = G V is stationary or equivariant, i.e. q(TjU),v) = q{uj^~ l v). An 

element w of the probability space will be called configuration. The mathematical 
expectation associated to the probability P will be denoted by E. 

For a configuration uj, a site w with q(ui, v) = will be called active or undeleted 
and a site t> with t>) = 00 deleted. 

For each W G fi denote by 

V(w) - V(X(w)) = {w G V I w) = 0} 

the subset of active vertices, and denote by A(w) the corresponding induced sub- 
graph of X. It is the sub-graph of A whose vertex set is V(u>) and whose edge set 
is 

E(u) = E(X(uj)) = {e G E(X) \ both cndpoints of e arc in V(w)} 

Let A = (y(A), E(A)) be an (deterministic) induced sub-graph of A. It gives 
rise to a random family of induced sub-graphs A(w) := A(w) n A. 

On any of the graphs introduced so far we will consider operators of finite hopping 
range. The easiest example to have in mind is the adjacency operator A u considered 



4 



IVAN VESELIC 



already in (1). More generally, a operator of finite hopping range H on a graph G 
is a linear map H : 1 2 {V(G)) — > l 2 (V(G)) such that there exists C, R < oo with 

(i) H(v,w) = H(w,v) 

(ii) H(^v,^w) — H{v,w) for all 7 E Y 
(hi) < C and 

(iv) H(v,w) = if dist(v,w) > R 

for all v,w E V(G). Here H(v,w) :— (8 V ,H5 W ) and <5„ E l 2 (V(G)) is the function 
taking the value 1 at v and elsewhere. 

For a sub-graph G C X and a finite hopping range operator H denote by H G 
the compression of H to l 2 (V(G)), in other words 

H G (v, w) = H(v, w) if v, w E G and H G (v, w) = otherwise 

If V = F(G) is finite, is a (|F| x |y|)-matrix, where | • | denotes the cardinality of 
a set. Thus the spectrum of H G is real and consists entirely of eigenvalues Xi(H G ), 
which we enumerate in increasing order counting multiplicity. Let us define the 
normalized eigenvalue counting function of H G as 

N (H,X).= — 

We assume that the discrete group T is amenable, i.e. there exists a F0lner 
sequence {Ij}j of finite, non-empty subsets of T. A sequence {Ij}j is called F0lner 
sequence if for any finite K C T and e > 

(3) \IjAKIjl < e\Ij\ 

for all j large enough. Since the quotient X/Y is compact, it follows that K :— {7 E 
r I jTDJ 7 7^ 0} is a finite generator set for T, cf. §3 in [2] for a similar statement in 
the context of manifolds. Now for finitely generated amenable groups there exists 
a F0lner sequence of subsets, which is increasing and exhausts T, cf. Theorem 4 in 
[1]. From [15] we infer that each F0lner sequence has an tempered subsequence. A 
tempered F0lner sequence is a sequence which satisfies in addition to (3) the growth 
condition 

there exists G < 00 such that for all j E N : < G\Ij\ 

To each increasing, tempered F0lner sequence associate an admissible exhaustion 
{Aj}j of X given by 

A,:- |J 7^CA 
lei; 1 

where IJ 1 := {7I7- 1 s Ij}. 

For a finite hopping range operator H, a F0lner sequence {Ij}j, and a random 
configuration uo E £1 introduce for brevity sake the following notation: :— 
H x ^\ HI := H A i("\ and N^(X) := N(H^X). Denote by P U (I) := xi(H^) the 
spectral projection of associated to the energy interval /. 

Theorem 2.1. There exist a distribution function N called integrated density of 
states such that for almost all u E ft and any admissible exhaustion Aj , j E N we 
have 



(4) 



lim Nl{E) = N(E) 

j—>oo 
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at all continuity points of N . The following trace formula holds for the IDS 
N(E) = ±E{Tr( X rP u (\-oo,E[))} 

We say that the IDS N is associated to the sequence of random operators 
{Hi}uen,j £ N. Next we address the question of boundary condition dependence. 
Denote A c = X \ A. 

Proposition 2.2. Let H be a finite hopping range operator, Aj,j G N an admissible 
exhaustion and R G N, C < oo. Let B 3 ' : l 2 (Aj) — ► l 2 (Aj),j G N be any sequence 
of symmetric operators such that for all v,w G V we have \B J (v,w)\ < C and 
B j (v,w) = i/dist(w,Ap +dist(io,Ap > R. Then the IDS' associated to the 
sequences {i?i} we fi, j G N and {Hi + B 3 } uen ,j G N coincide. 

Next we establish the non-randomness of the spectrum of H u and its compo- 
nents, its relation to the IDS and an understanding of the IDS as a von Neumann 
trace. Denote by a disc, cr ess , a ac , c sc , a pp the discrete, essential, absolutely continu- 
ous, singular continuous, and pure point part of the spectrum. Denote by cr comp the 
set of eigenvalues which posses an eigenfunction with compact, i.e. finite, support. 
In the following theorem T need not be amenable, but X must be countable. 

Theorem 2.3. There exists a O' C O of full measure and subsets of the real 
numbers S and S #; where • G {disc, ess, ac, sc,pp, comp}, such that for all wed' 

a(H u ) = H and <r.(H u ) = -£ m 

for any • = disc, ess, ac, sc,pp, comp. IfT is infinite, S^isc = 0- The almost-sure 
spectrum S coincides with the set of points of increase of the IDS 

S = {A G R | N{\ + e) > N(X - e) for all e > 0} 

Furthermore, N is the distribution function of the spectral measure of the direct 
integral operator 

H:= H^aTiuj) 
Jn 

On the von Neumann algebra associated to Ti. there is a canonical trace and N(E) 
is the value of this trace on the spectral projection of H associated to the interval 
]-oo,E[. 

3. Proofs of the theorems 

Let H be a finite hopping range operator and assume without loss of generality 
\H(v, w)\ < 1 for all matrix elements. It follows that the / 2 -norm of H is bounded 
by K := 2dcg^. Since H is symmetric, it is a selfadjoint operator. In particular 
the spectrum of H w is contained in [-K, K] for all u) G O. 

Each 7 G r induces an unitary operator £7 Wj7 : ^{V^^uj)) — > 1 2 (V(uj)), (C/ w , 7 /)(«) := 
f(-f~ 1 v). Note that V^uo) = -fV(iv). By the definition of t 7 the action of T on ft 
and on X is compatible: 

(5) U^.-yH^U^^ — H TiUJ 
The equivariance formula (5) implies 

(6) U u „f(H u )UZ„ - f(H TyU ) 
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for any polynomial /. For continuous functions /, g we have \\f(H u ) — q(H u! )\\ < 
||/ - ffHoo. ' Thus /„ - / in || • ||oc) implies f n (H u ) -> in op- 

erator norm, and (6) extends by Weierstrafi' approximation theorem to all / E 
C([—K,K]). By taking scalar products we obtain the corresponding equivariance 
relation for the matrix elements: 

/(^)( 7 - 1 U ,7- 1 W ) = f(H Ty0J )(v,w) 



For the proof of the main Theorem 2.1 we need two key ingredients: an estimate 
of boundary effects on traces and a sufficiently general ergodic theorem, which will 
be applied to trace functional of the type 

F(w) := J2 W u ){v,v) = l^" 1 Tr (f(H u ) X r) 

Let us first estimate the boundary effects. 
Proposition 3.1. Let f(x) = x m for m E N. Then 

sup \Tr(f(Hl)) - Tr( XA ,/(^))| - 
wen |Aj| 

as j — > oo. 

Proof. We introduce the notion of a thickened boundary on a graph. For a sub- 
graph A and h E N set c^A := {v E A dist(w, A c ) < /i}. We expand the trace of 
powers of H^: 

Tr(Hir=J2( H ti m (v,v)=Y, E H u {v,v 1 )...H u {v m - 1 ,v) 

By an analogous formula for Tr{xAjH™) we obtain 

Tr[ XAj H™-(Hir]= J2 '£ l H u (v,v 1 )...H u (v m - 1 ,v) 
veAj • 

< I^A^degf* 

where the bullet denotes summation over m — 1-tuples (paths) in V(X) with at 
least one vertex outside Aj . By the F0lncr property of the sequence Ij , j E N we 
have: 

I^A,| 



lim f = for any h > 

j^OC \Aj\ 

This is the content of Lemma 2.4 in [17]. In fact, there manifolds are considered, 
but the proof applies literally to the case of graphs. □ 

Lindcnstrauss proved in [15] a remarkable ergodic theorem which applies to lo- 
cally compact, second countable, amenable groups. It includes the following state- 
ment for discrete groups. 

Theorem 3.2. Let T be an amenable discrete group and (CI,Bq,¥) be a probability 
space. Assume that T acts ergodically on by measure preserving transformations 
r 7 . Let {Lj}j be a tempered F0lner sequence in T. Then for every F E 



(7) lim -L £ F(t 7 w) = E{F} 

for almost all w E il. 



3 1 -yeij 
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In the application we have in mind F £ L°°, so the convergence holds in the 
i 1 -topology, too. 



Proof of Theorem 2.1. To prove the claimed convergence of distribution functions 
it is sufficient to establish f fdN£ -> / fdN for all / e C([—K, K]). This can in 
turn be reduced to proving the convergence in the case where / is a polynomial: 
By Weierstrafi' approximation theorem the polynomials are dense in C{\—K,K\). 
Let fk be a sequence of polynomials such that \\fk — f\\oo —* 0. For any e > one 
can choose k large enough such that ||/fc — /Hoc < e/3 and subsequently j large 
enough (depending on e and fk) such that \J fkdN£ — J fkdN\ < e/3. Since the 
total measure of N and any N£ is bounded by one, 



/«-/. 



fdNi - I fdN 



< 



+ 



f fkdNZ - J f k dN + J (f k - f)dN 



fk)dNi 

Thus it is sufficient to prove the convergence of moments 

(8) lim / \ m Nl{d\)= ( X m N(dX) for all m e 



< e 



N 



Next we show that the limit on the left hand side equals -^E {Ti(xrH" 1 )} f° r 
almost all u. For this aim we write the moment of the IDS as a trace 



/ \ m Ni{d\) = \K j \- 1 Tr{f{Hi)) 



which by Proposition 3.1 converges for j — > exo to the same limit as |Aj| _1 Tt(xa., ,f{H w )). 
Now we decompose the trace according to local contributions and apply Lindcn- 
strauss' theorem 

lA^TrixAjiH^)) = \Aj\~ 1 ]T f(H> u ){v,v) = 

m~ x E E m^-'v^-'v) = \ij\- 1 J2 i^r 1 E 

— > E{F} as j — > oo for almost all w e H 

where F(w) = l^l" 1 £^ f(H„)(v,v) = l^" 1 Tr( X ^/(if w )). Set 
Eu,(A) := P w (] — oo, A[). The expectation of F equals 



1 1 l^e^ 



A m F w (dA) («,?;) 



^ | A m E{Tr( X ^(dA))} = | A m iY(dA) 



□ 
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Proof of Proposition 2.2. Denote := + BK Similarly as in the proof of 
Proposition 3.1 we have 

Tr[(^r - mn 

= ^ ^H u {v,v{) . . .H u (v m -i,v) - H u (v,vx) . . .H^Vm-i^) 

<{l + (l + Cr]\d Rm+it A 3 \degfr 

Here the bullet denotes the summation over all paths v± , . . . , v m -i E Aj with at least 
one vertex in d^Aj and r := max(i?, R). By the F0lner property |Aj| _1 Tr[(i^) m — 
(H^) m ] converges to zero as j ^> oo. □ 

Proof of Theorem 2.3. First we prove the non-randomness of o~ com p- Set 
(9) £ :={E G M | 3 finite induced sub-graph G C X and / G ^ 2 (G) 

such that i? G / = £/} 

Then o- CO mp{H-u) C X for all u; G fi. Since X is countable, there exists a countable 
exhaustion of X by finite sets Dj,j G N. If we set 

:={£:el|3weO and / G J 2 (Dj(a/)) s.t. H D ^ f = Ef}, 

then E = UjSj and thus E is countable. 

For any £ e S set Ob := {cj | 3/ with finite support and H^f — Ef}. This 
set is invariant under the ergodic action of T by the transformations r 7 . Therefore, 
either P(Ob) = 1 or P(Q.e) = 0. In the first case set £Ie = He and -E G E comp , 
in the second set He = H E an d E & ^comp- Here the superscript c denotes the 
complement of a set. The set f2 := dggj, ^ has full measure and each iv E ft 
satisfies o~ corri p[B Lt j} E com p. 

The remaining statements of the theorem follow from the results of [13]. One 
just has to check that the required assumptions are satisfied. This is not hard, but 
it would require to introduce the notion of grupoids and basic features of Connes' 
non-commutative integration theory [6] . Therefore we leave the details of the proof 
of Theorem 2.3 to another occasion. □ 

4. Outlook: finitely supported and exponentially decaying states 

Once we have a rigorous definition of the integrated density of states for the 
QPM, we can study finer spectral properties. One of the main interest in the 
physics literature are the properties of bound states, and their contribution to nu- 
merically observed "peaks" of the density of states. This quantity is the distribu- 
tional derivative of the IDS. In the following we restrict our discussion to the QPM 
corresponding to the adjacency operator A on the lattice Z d , and to the sequence 
q(-,v),v E Z d , consisting of independent, identically distributed random variables. 

There seem to be three different types of bound states of the QPM: finite cluster 
states, molecular states and exponentially decaying states, cf. [8, 7, 12, 18, 5]. Finite 
cluster states occur since almost surely there are active clusters of finite size, which, 
consequently, can support only bound states. These states are mathematically not 
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challenging. 

However, as pointed out earlier, there exist molecular states. They are eigen- 
vectors of the adjacency operator restricted to the active sites with support on a 
finite region of the infinite cluster. This means that due to the deletion of sites the 
unique continuation property of eigenfunctions of the adjacency operator on the 
lattice breaks down. 

For the analysis of molecular states it is convenient to introduce the restriction 
of A to the infinite active cluster X°°(u), and to define the corresponding IDS 

by 

N°°(\) = |J r | _1 E{Tr[x_ ? roo( w )£'2 (A)]} 

where T°°(uj) = T^X°°{uj) and T = {0}. N°° is self-averaging, i.e. can be defined 
by an exhaustion procedure, similarly as N in Theorem 2.1. 

Here is a result on molecular states, whose proof is given in [20]. 

Theorem 4.1. The set of points of discontinuity of the IDS N of {A^}^ coincides 
with the set of points of discontinuity of the IDS N°° of {A^}^ and equals S as 
defined in (9). 

Finally, one can hope that a multi-scale argument will yield a proof of the exis- 
tence of exponentially decaying eigenstates, as in the case of the Anderson model 
[3, 9, 21, 19]. For this aim, one has to develop new tools to deal with the singular 
randomness present in the QPM. On the other hand, once this is done, one might 
use the new ideas to approach the exponential localization problem for the Anderson 
model with Bernoulli disorder of the coupling constants. In the multi-dimensional 
case this is a problem which is open since decades. 
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